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ABSTRACT 

Measurements of neutron star masses and radii are instrumental for deter- 
mining the equation of state of their interiors, understanding the dividing line 
between neutron stars and black holes, and for obtaining accurate statistics of 
source populations in the Galaxy. We report here on the measurement of the 
mass and radius of the neutron star in the low-mass X-ray binary KS 1731—260. 
The analysis of the spectroscopic data on multiple thermonuclear bursts yields 
well-constrained values for the apparent angular area and the Eddington flux 
of the source, both of which depend in a distinct way on the mass and radius 
of the neutron star. The binary KS 1731—260 is in the direction of the Galactic 
bulge, allowing a distance estimate based on the density of stars in that direction. 
Making use of the Han & Gould model, we determine the probability distribu- 
tion over the distance to the source, which is peaked at 8 kpc. Combining these 
measurements, we place a strong upper bound on the radius of the neutron star, 
R < 12 km, while confining its mass to M < 1.8 M . 

Subject headings: stars: neutron — X-rays: binaries — stars: individual (KS 1731— 260) 



1. Introduction 

A direct probe of the equation of state of cold, ultradense matter is the measurement 
of the radii of neutron stars, whose cores are comprised of the densest matter in the current 
universe. Radii measurements can distinguish between a variety of interior compositions, as 
well as the strength of the many-body nuclear force (Lindblom 1992; Lattimer & Prakash 
2001; Read et al. 2009; Ozel & Psaltis 2009). In addition, the mass distribution of neutron 
stars can give insights to the outcomes of supernova explosions and evolutionary tracks in 
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binaries. If sufficiently high, a mass measurement can also lead by itself to constraints on the 
composition and the interactions of the interior (Demorest et al. 2010; Ozel et al. 2010a). 

While precise mass measurements of neutron stars have been possible for decades, sig- 
nificant progress in the determination of radii of neutron stars has been made only recently. 
There are, by now, a handful of neutron stars for which both a radius and a mass measure- 
ment have been achieved. Some of these spectroscopic observations have been carried out 
during the quiescent episodes of accreting sources that reside in globular clusters (Rutledge 
et al. 2001; Heinke et al. 2006; Webb & Barret 2007; Guillot et al. 2010). They have led to 
significant, albeit correlated, constraints in the mass and radius of neutron stars. 

Another way to measure the neutron star mass and radius is through a combination of 
spectroscopic phenomena observed from their surfaces during thermonuclear X-ray bursts 
(van Paradijs 1979; Ozel 2006). The repeated high count-rate bursts allow a measure of 
the apparent angular area of the neutron star over a wide range of temperatures using 
time-resolved spectra. Furthermore, very luminous X-ray bursts (called photospheric radius 
expansion, or PRE, bursts) cross the local Eddington flux at the neutron star surface and 
allow us to obtain a measure of the neutron star mass (as corrected for general relativistic 
effects). Combined with the distance measurement to the source, these quantities can be 
converted into independent measurements of the neutron star mass and radius. This ap- 
proach led to the determination of the masses and radii of several neutron stars (Ozel et al. 
2009; Giiver et al. 2010a, b) and enabled the measurement of the pressure of neutron star 
matter above nuclear saturation density (Ozel et al. 2010b). 

In this paper, we place strong constraints on the mass and radius of KS 1731—260 based 
on the analysis of thermonuclear bursts observed with the Rossi X-ray Timing Explorer. 
The transient X-ray burster KS 1731— 260 was discovered in 1989 during observations with 
the TTM/Kvant telescope on board the Mir station (Sunyaev et al. 1990). Three Type 
I bursts were seen during the immediate follow-up observations, confirming the presence 
of a neutron-star in this X-ray binary. RXTE detected a total of 24 Type-I bursts from 
KS 1731—260, which allow a systematic study of the spectra obtained during the bursts and 
the measurement of the surface area and the Eddington limit of this neutron star (Galloway 
et al. 2008; Giiver et al. 2010c, d). 

In order to convert the measurement of the apparent angular area and Eddington flux 
from the surface of a neutron star into a measurement of its mass and radius, we need an 
estimate of the distance to the source. In earlier work, we used sources in globular clusters 
with known distances (Ozel et al. 2009; Giiver et al. 2010b) or measured the source distance 
using the method of red clump stars (Giiver et al. 2010a). In this paper, we use the fact 
that KS 1731— 260 lies within the Galactic bulge, which has a limited spatial extent in the 
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Galaxy. Even though the further localization of the source within the bulge is not easy to 
achieve, it nevertheless allows us to place a strong upper bound on the mass and radius of 
the neutron star. 

In Section 2, we present the constraints on the distance to KS 1731—260. In Section 3, 
we summarize the results from the analysis of the spectra of thermonuclear bursts from this 
source and show the resulting constraints on its mass and radius. In Section 4, we compare 
these measurements to other existing mass-radius measurements. 



2. The Distance to KS 1731-260 



, i.e., almost exactly at the 



KS 1731— 260 lies at galactic coordinates (7,6) = (1.06,3.65) 
position of Baade's window (reflected in the Galactic plane). Its position toward the bulge 
allows an estimate of its distance based on the distribution of stars in that direction. Having 
no prior information on the distance to KS 1731—260, we will assume that the likelihood 
that it resides at a given distance D is proportional to the number density of stars along the 
line of sight to the bulge at that distance. 



For the stellar distribution along the line of sight we adopt the model of lHan fc Gould 



( 120031 ). which is based primarily on star counts, and, without any adjustment, reproduces the 



microlensing optical depth measured toward Baade's window. For the b ulge. lHan fc Gould 



( 120031 ) normalize the "G2" K-band integrated-li ght-based bar model of iDwek et al.l (11 995) 



using star counts toward Baade's window from iHoltzman et al.l (119981) and IZoccali et al. 
( 120001 ). For the disk, they incorporate the model of IZheng et al.l (120011 ). which is a fit to 
star counts. In this model, there is a 1 kpc hole at the center of the Galactic disk, since 
this material would have gone into the Galactic bulge according to the standard picture of 
pseudo-bulge forma tion. We normali ze the distance scale to a Galactocentric distance of 
R = 8.0 ± 0.4 kpc (jYelda et al. 2010 ). The position of the source, being almost exactly at 
the position of Baade's window (reflected in the Galactic plane), minimizes the extrapolation 
from calibrating star counts to the adopted model. 

In Figure 1, we plot the stellar distribution along the line of sight to KS 1731—260, which 
we then take to be equal to the likelihood of the source distance. The dotted lines show 
the disk and bulge contributions separately, while the solid line shows the total likelihood 
for a distance to the Galactic center of 8.0 kpc. The dashed line shows the likelihood of the 
distance to KS 1731—260 when a Gaussian uncertainty of 0.4 kpc has been incorporated to 
the distance to the Galactic center. It is evident from the dashed curve that the details of the 
prescription for the 1 kpc hole in the Galactic disk are not important. The latter likelihood, 
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Fig. 1. — The solid line shows the likelihood of the distance to KS 1731—260, which we have 
taken to be proportional to the number density of stars along the line of sight to the source, 
at each given distance, for an assumed distance to the Galactic center of 8.0 kpc. The dotted 
lines show the disk and bulge components of this distribution. The dashed line shows the 
total likelihood of the distance to KS 1731—260 when a Gaussian uncertainty of 0.4 kpc has 
been incorporated to the distance to the Galactic center. 
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which we use hereafter, indicates that the source most probably lies at 5 kpc< D <11 kpc. 



3. Determination of the Neutron Star Mass and Radius from X-ray Burst 



RXTE observed KS 1731-260 for a total of 483 ks until June 2007. We reported in Giiver 
et al. (2010c) the systematic analysis of the 24 bursts detected in this time period. All of the 
details related to the extraction of the spectra, detector response, background subtraction, 
and the spectral fitting are discussed in Giiver et al (2010c). We only summarize here the 
results of the analysis of the time-resolved continuum spectra in order to utilize them for the 
mass and radius determination of the neutron star in this binary. 

The 24 bursts we use here met our requirement that the persistent flux prior to the burst 
was < 10% of the peak burst flux, which minimizes uncertainties related to background 
subtraction. We analyzed a total of 1309 spectra from these bursts, out of which 1240 
spectra (98%) gave acceptable values of reduced x 2 . We measured the evolution of the flux, 
the temperature, and the angular size using the spectra obtained from 0.25, 0.5, or 1 s 
integrations during each burst. The bursts follow highly reproducible tracks on the flux- 
temperature diagram, as is shown in Figure 2. The slight decline in the normalization at 
very low fluxes is consistent with the expected increase of the color correction factor at low 
temperatures. Incorporating this trend as a systematic uncertainty, Giiver et al. (2010c) 
reported an apparent angular area of A = 88.4 ± 5.1 (km/10 kpc) 2 , which corresponds to a 
radius of 9.4 ± 0.3 (km/10 kpc) during the cooling tails of bursts. The quoted uncertainty 
in the radius includes statistical and systematic uncertainties. The apparent angular area 
depends on the stellar parameters according to 



where f c is the color correction factor. 

Two out of the 24 bursts from KS 1731—260 showed clear evidence of photospheric radius 
expansion, where the photosphere in highly energetic burst events expands to many times the 
stellar radius, while the flux and temperature follow a characteristic evolution (see Giiver et 
al. 2010d). We measured the Eddington flux in these events at the touchdown moment when 
the photosphere returns to the stellar surface, as indicated by a maximum in the temperature 
and a minimum in the inferred area (as in Ozel et al. 2009 and Giiver et al. 2010d). For 
KS 1731— 260, we found an average touchdown flux of (4.45 ± 0.12) x 10~ 8 erg cm" 2 s -1 , 
where the quoted uncertainty is purely formal. Because in sources with larger number of 
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Fig. 2. — The evolution of the bolometric flux and the blackbody temperature during cooling 
tails of three thermonuclear X-ray bursts observed from KS 1731—260. The bursts follow the 
L oc T 4 relation, which is indicated by the solid line. 
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Fig. 3. — The confidence contours over the blackbody temperature and normalization during 
the touchdown moments of the two photospheric radius expansion bursts (in green and 
black, respectively) observed from KS 1731—260. The two bursts differ in both temperature 
and angular area, and they each have significant uncertainties in these variables, which are 
strongly correlated. However, when the distribution is projected onto Eddington flux (red 
curves), the result is quite consistent between the two bursts and well constrained. 
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photospheric radius expansion bursts we found systematic uncertainties of the order of ~ 5% 
in the touchdown flux, we assign hereafter a total uncertainty of 0.22 x 10~ 8 erg cm -2 s _1 
in order to be conservative. The touchdown flux depends on the stellar parameters as 



where G is the gravitational constant, c is the speed of light, and k es is the opacity to electron 
scattering. 

We assign independent probability distribution functions to the distance P(D)dD, the 
touchdown flux P(F<td)^-Ftd, and apparent angular area normalization P(A)dA based on 
the measurements and the distance estimate discussed above. We take the probability distri- 
butions over the touchdown flux and the apparent angular area to be Gaussian. Converting 
the Eddington limit and the apparent angular area to a mass-radius measurement requires 
a prior knowledge of the hydrogen mass fraction X in the atmosphere and a model for the 
color correction factor f c . In the absence of any information on the composition of the ac- 
creted material for KS 1731—260, we take a box-car distribution for X that covers the range 
0—0.7. The color correction factor, on the other hand, is obtained from modeling the hot 
atmospheres of accreting, bursting neutron stars (e.g., Madej, Joss, & Rozanska 2004). We 
take the box-car probability distribution centered at f c0 = 1.35 and a width Af c = 0.05 that 
is appropriate for a thermal flux in the range between « 1% — 50%, which is seen in the 
cooling tails of bursts (see Giiver et al. 2010b, 2010c for a detailed discussion). 

We calculate the probability distribution over the neutron star mass and radius as in 
Giiver et al. (2010a). Note that we have corrected here an error in the Equation (5) of 
Ozel et al. (2009) where the term 7GM/Rc 2 should read 8GM/Rc 2 . Figure 1 (left panel) 
shows the 68% and 95% confidence levels for the mass and the radius of the neutron star 
in KS 1731—260. Although the observed apparent angular area and the touchdown flux are 
consistent with a range of masses and radii for the neutron star, they nevertheless provide 
a strong upper bound of R < 12 km and M < 1.8M Q . In Figure HI we also identify the 
astrophysically relevant range of neutron star masses M > 1.2M and show the combined 
constraints as filled contours. The allowed range of neutron star radii is inconsistent with 
equations of state that predict large neutron star radii, ~ 12 — 15 km. This is in agreement 
with the earlier spectroscopic measurements of radii in bursting neutron stars (Ozel et al. 
2009; Giiver et al. 2010a, b) as well as the quiescent neutron stars in globular clusters (Webb 
& Barret 2007; Guillot et al. 2010). 

When performing parameter estimation in models within a Bayesian framework, two 
concerns often arise. First, there is no universally accepted measure of goodness-of-fit, which 
makes it difficult to assess whether data sets that are used to estimate the parameters are 
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Fig. 4. — (Left) The 68% and 95% confidence contours for the mass and radius of the neutron 
star in KS 1731— 260. The astrophysically relevant range of masses M > 1.2M are shown 
as filled contours. The representative mass-radius relations for a select number of equations 
of state are also shown and include multi-nucleonic ones (A, FPS, UU, AP3), equations of 
state with condensates (GS1-2), strange stars (SQM1, SQM3), and meson-exchange models 
(MSO). The black line indicates the black hole event horizon. The descriptions of the various 
equations of state and the corresponding labels can be found in Lattimer & Prakash (2001) 
and Cook, Shapiro & Teukolsky (1994). (Right) The 68% and 95% confidence levels over 
neutron-star mass and radius calculated by using the maximum combined likelihood along 
lines of constant distance and hydrogen mass fraction rather than by integrating over these 
two nuisance parameters. The fact that these confidence levels are very similar to those 
shown in the left panel demonstrates that marginalization does not bias the results. 
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consistent with the model. Second, in models with more than two parameters, marginalizing 
(integrating) over the so-called nuisance parameters in Bayesian analysis in order to infer the 
parameters of interest can introduce biases. This is especially a concern when measurement 
uncertainties in different parameters differ widely from each other. In our case, the allowed 
logarithmic range of values over distance is significantly larger than the uncertainties over 
the touchdown flux and the apparent radius. These concerns have led Steiner et al. (2010) 
to explore a different, ad hoc, interpretation of the Eddington flux in their inference of the 
neutron star masses and radii. 

In Figure [5j we show the contours in the mass-radius plane that correspond to each 
individual spectroscopic measurement, taking into account the la ranges in the measured 
quantities but fixing the distance to a representative value of 7 kpc, which is within the region 
of highest likelihood (see Fig. [Q and in the range of distances most favored by the data. The 
two sets of contours intersect with each other within their la levels and overlap with the 
combined confidence contours over mass and radius. Even though this is not a quantitative 
measure of goodness-of-fit, it shows the self-consistency of the theoretical framework with 
the data. 

In order to address the second concern about the effect of marginalizing over distance, we 
performed the following test. We calculated the likelihood over mass and radius by using the 
maximum of the combined likelihood along lines of constant distance and constant hydrogen 
mass fraction X, instead of integrating over these two parameters. We show in the right 
panel of Figure H] the resulting confidence contours over mass and radius. The result is very 
similar to that shown in the left panel of Figure HJ demonstrating that the marginalization 
over the two most uncertain parameters, namely distance and the hydrogen mass fraction, 
does not bias the mass and radius measurement. 



4. Conclusions 

We used the spectroscopic measurements of the apparent angular area and the Ed- 
dington flux during thermonuclear X-ray bursts from KS 1731—260 to measure its mass and 
radius. We place a strong upper bound on the radius of R < 12 km and confine the mass to 
M < 1.8M . This measurement challenges equations of state that predict radii larger than 
~12 km. 

We also explored whether, within the Bayesian parameter estimation approach, 1. the 
theoretical framework and the data sets used here are self-consistent and 2. the integration 
over the nuisance parameters introduces biases in the confidence limits over the neutron star 
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Fig. 5. — The contours in the mass-radius plane that correspond to the measurement of the 
Eddington flux and the apparent radius for a representative distance of 7 kpc. The dashed 
lines reflect the la errors on the measured quantities and model parameters. The fact that 
the two sets of contours intersect at the location of the confidence limits shown in the left 
panel of Figure @] and overplotted here (filled contours) demonstrate the consistency of the 
model with the data. 
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mass and radius. We found that for a highly likely value of the distance, the contours in 
the mass-radius plane that correspond to the apparent radius and the Eddington limit are 
consistent with each other and with the confidence contours over mass and radius within 
la. Moreover, we showed that considering the most likely configuration for each value 
of the distance and the hydrogen mass fraction, as opposed to integrating over these two 
parameters, leaves the results largely unaffected. These results justify our identification of 
the touchdown point for the measurement of the Eddington flux and our conclusion that the 
apparent angular area during the cooling tail of the burst corresponds to the entire neutron 
star surface. 

This work was supported in part by NASA ADAP grant NNX10AE89G and Chandra 
Theory grant TMO-11003X. AG was supported by NSF grant AST-0757888. We thank 
Andrew Steiner for pointing out an error in a numerical factor in our earlier calculation 
of the Jacobian transformation. We thank Dimitrios Psaltis for a careful reading of the 
manuscript. 
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